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Abstract 

It has been shown that the evolution of collectivity in ultra-relativistic heavy ion 
collision is manifested in the variation of various HBT radii with invariant mass 
(M) extracted from the correlation functions of two lepton pairs. The value of 
the radial velocity (v r ) can be estimated from the ratio of the px distributions 
of single photons to lepton pairs for various M windows. It has been argued that 
the variation of radial flow with appropriate kinematic variables can be used as 
an indicator of a phase transition from initially produced partons to hadrons. We 
also consider the elliptic flow {v§ F ) of the matter as probed by the single electron 
spectra originating from the semileptonic decays of heavy mesons. The measured 
values of v^ F and the nuclear suppression factor (Raa) a t RHIC energy have been 
reproduced simultaneously by including both the collisional and radiative processes 
within the scope of perturbative quantum chromodynamics. The Raa and v^ F have 
been predicted for LHC energy. 
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1 Introduction 



The hot and dense matter expected to be formed in the partonic phase after ultra- 
relativistic heavy ion collisions (URHIC) dynamically evolve in space and time due 
to high internal pressure. The system cools due to expansion and reverts back to 
hadronic matter from the partonic phase. It is well known that the average mag- 
nitude of radial flow at the freeze-out surface can be extracted from the transverse 
momentum (py) spectra of the hadrons. However, hadrons being strongly interacting 
objects can bring the information of the state of the system when it is too dilute to 
support collectivity. The electromagnetic (EM) probes, i.e. photons and dileptons 
on the other hand are produced and emitted p] (see [2] for review) from each space 
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time points. Therefore, estimating radial flow from the EM probes will shed light on 
the time evolution of the collectivity in the system. The generation of collectivity in 
the system depends on EoS - hence this can be used to differentiate partonic and 
hadronic phases as the EoS for these two phases are different. In case of EM probes- 
dilepton has the advantage over the real photons. Because the low pt photons from 
the hadronic phase receive large transverse kick due to radial flow and consequently 
appear in the high pt domain to mingle with those from the QGP phase, making the 
detection of photons from QGP difficult. However, for dileptons there are two kine- 
matic variables available - the pt and the invariant mass (M). While the pt spectra 
of dilepton is affected by the flow, the pt integrated M spectra remains unaltered. 
This suggests that a careful selection of px and M windows will be very useful to 
characterize the QGP and hadronic phases. In the present work we will demonstrate 
how the development of radial flow can be estimated through the Hanbury-Brown 
Twiss (HBT) interferometry with virtual photons (lepton pairs) for different M win- 
dows. The radial flow velocity (v r ) can be estimated by considering the ratio of the pt 
distribution of single photon to lepton pairs for various M windows. We will briefly 
discuss the results here and refer to [3] for details. 

Single electrons originating from the semi-leptonic decays of heavy mesons carry 
the information on the interaction of the heavy quarks (a constituent of the heavy 
mesons) with the thermal medium of light quarks and gluons produced in heavy ion 
collisions. The Raa and vf F can be used to quantify the interaction of the heavy 
quarks with the QGP. Several ingredients like inclusions of non-perturbative contri- 
butions from the quasi-hadronic bound state [3], 3-body scattering effects [5], the 
dissociation of heavy mesons due to its interaction with the partons in the ther- 
mal medium [6] and employment of running coupling constants and realistic Debye 
mass [7], the inclusion of both elastic and inelastic collisions along with the path 
length fluctuation have been proposed [8J to improve the description of the experi- 
mental data. Within the framework of Fokker Planck equation (FPE) we will evaluate 
v^ F and Raa for these electrons. In the present paper we discuss the elliptic flow 
of the matter probed by the single electron from the heavy mesons decays. For the 
elliptic flow of the matter probed by single photon and lepton pair we refer to PfTU] 
for details. 

In the next section we will briefly describe the HBT interferometry with virtual pho- 
tons. The azimuthal anisotropy of the system probed by single electrons originating 
from the heavy flavour decays will be discussed in section 3. Section 4 is devoted to 
summary and discussions. 



2 HBT interferometry with dileptons 

The interferometry of the dilepton pairs actually reflect correlation between two 
virtual photons, the analysis then can proceed by computing the Bose-Einstein cor- 
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Fig. 1. Left panel: variation of R s ide and Rout as a function of (M). The dashed (solidi 
with asterisk) line indicates HBT radii for QGP (total=QGP+hadron) phase. Right panel: 
variation of the radial velocity with (M) extracted from the ratio of the pt distribution of 
photons to lepton pairs (see [3] for details). 

relation (BEC) function for two virtual photons which can be denned as, 6*2(^1,^2) = 
P2(pl,P2)/ [Pi(pi)Pi{p2)), where pi is momentum of the individual lepton pair, P\{pi) 
and P2{p~i,P~2) represent the one- and two- particle inclusive lepton pair spectra re- 
spectively, which can be evaluated form the source function for various invariant mass 
windows of the pair [11]. For the productions of lepton pairs from QGP the anni- 
hilation of thermal quarks and from the hadronic phase the decays of thermal light 
vector mesons (p, u and 0) have been considered. 

For the space time evolution of the system relativistic hydrodynamical model with 
cylindrical symmetry [12] and boost invariance along the longitudinal direction [13] 
has been used. For a system undergoing isentropic expansion, the initial temperature 
(Tj) and proper thermalization time (rj) of the system may be constrained by the 
measured hadronic multiplicity, dN/dy ~ Tfr^. For Relativistic Heavy Ion Collider 
(RHIC) we have taken Tj = 290 MeV and Tj = 0.6 fm/c. The EoS which controls the 
rate of expansion/cooling has been taken from the lattice QCD calculations [14]. 
The chemical (T c ^=170 MeV) and kinetic (T^=120 MeV) freeze-out temperatures 
are fixed by the particle ratios and the slope of the pt spectra of hadrons [T3] . With 
all these ingredients the correlation function C2 has been evaluated for different 
(average) invariant mass windows, (M) = (Mi + M^jl as a function of q S ide and 
q ou t [11] which are related to the transverse momentum of individual pair. The HBT 
radii, R S id e and R out corresponding to q si d e and q out extracted from the (Gaussian) 
parametrization of Ci. 

The R S ide is related to the transverse size of the system whereas the R out mea- 
sures both the transverse size and duration of particle emission ( [TB] for review). 
The Rside shows non-monotonic dependence on (M) (Fig. (TJ left panel). It can be 
shown that R side ~ 1/(1 + E collective / E therma i). The high (M) regions are dominated 
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by the early partonic phase [16] where the collective flow has not been developed 
fully consequently the ratio of collective (E coUective ) to thermal (E thermal ) energies 
is small- hence a larger R S id e is obtained for large M. In contrast, the lepton pairs 
with M ~ m p are emitted from the late hadronic phase where the collective flow or 
1 + E coUective /E therma i) is large, which is reflected as a dip in R side for (M) ~ m p . 
Thus the variation of Rside with M can be used as an efficient tool to measure the 
collectivity in various phases of matter. We observe that by keeping the p and u 
contributions and setting radial velocity, v r — 0, the dip in R S id e vanishes, confirming 
the fact that the dip is caused by the large radial flow of the hadronic matter. The 
R ou t probes both the transverse dimension as well as the duration of emission and 
unlike Rside, Rout does not remain constant even in the absence of radial flow. The 
large (M) regions are populated by lepton pairs from early partonic phase where the 
effect of flow is small and the duration of emission is also small - resulting in smaller 
values of R ou t- For lepton pair from (M) ~ m p region the flow is large which could 
have resulted in a dip as in R S id e in this M region. However, R out probes the duration 
of emission too, which is large for hadronic phase. The larger duration overwhelms 
the reduction of R out due to flow in the hadronic phase resulting in a bump in R out 
in this region of (M). 

As mentioned before the v r can be estimated from the ratio of the p? spectra of real 
photons to lepton pairs. Fig. [T] (right panel) shows the variation of v r with (M) both 
for SPS and RHIC conditions. The individual spectra of photons and lepton pairs 
are constrained by the available experimental data [3]. The v r increases with M up 
to M — M p then drops. From the invariant mass spectra it is known that the low M 
(below p mass) and high M (above <fi peak) pairs originate from a partonic source [3] . 
The collectivity (or flow) does not develop fully in the QGP resulting in smaller values 
of v r at both low and high M regions. Lepton pairs for M ~ m p originate from the 
late hadronic source which are largely affected by the flow resulting in higher values 
of v r . In summary, the value of v r for M below and above the p-peak is small but 
around the p peak is large - the resulting behaviour is displayed in Fig. [1] (right 
panel). Similar non-monotonic variation of the effective slope parameter of the pt 
distribution of lepton pairs for various M windows is observed in [TTJ . The evolution 
of v r as observed in Fig. [1] (right panel) is responsible for such behaviour. 



3 Elliptic flow probed by single electron spectra 

The heavy flavors, namely, charm and bottom quarks, play a crucial role in charac- 
terizing the QGP (see also [H]). As the relaxation time is larger for heavy quarks 
than light partons, the light quarks and the gluons thermalize faster. Therefore, the 
propagation of heavy quarks through QGP may be treated as the interactions be- 
tween equilibrium and non-equilibrium degrees of freedom and the FPE provides an 
appropriate framework [TH] for such studies. In this work we would like to evaluate 
v^ F and Raa of heavy flavours within the framework of FPE and contrast the results 
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with the available experimental data. The evolution of heavy quarks momentum dis- 
tribution function, while propagating through the QGP are assumed to be governed 
by the FPE, which reads, 



dl = d_ 

dt dpi 



Mp)f + ^-[B l3 {p)f} 



(1) 



where the kernels and B 
for | p |— > 0, Ai — > -ypi and 
co-efficients respectively 



/ d 3 ku(p, k)ki B^ = J d 3 ku(p, k)kik 
— > DSij, where 7 and D stand for drag and diffusion 



%3 are given by, Ai 



The basic inputs required for solving the FP equation are the dissipative co-efficients 
and initial momentum distributions of the heavy quarks. The (effective) drag and 
diffusion coefficients have been evaluated by taking in to account both the collisional 
and radiative processes [20]. In the radiative process the dead cone and LPM effects 
are included. In evaluating the drag co-efficient we have used temperature dependent 
strong coupling, a s (T) [21J. The Debye mass, ~ g(T)T also a temperature dependent 
quantity used as a cut-off to shield the infrared divergences arising due to the ex- 
change of massless gluons. The initial momentum distribution of heavy quarks in pp 
collisions have been taken from the NLO MNR [22] code. The solution of the FPE 
for the heavy quarks is convoluted with the fragmentation functions of the heavy 
quarks to obtain the pt distribution of the D and B mesons. For heavy-quark frag- 
mentation function, the Peterson function has been used. The solution of the FPE 
has been used to predict the pr spectra of the D mesons by following the procedure 
similar to blast wave method [23], the result is compared with experimental data [23] 
(Fig. [2J left panel) which indicate that the present data can not distinguish between 
the equilibrium and non-equilibrium scenario. The px distribution of the electrons 
from the semi-leptonic decays of D and B meson are evaluated using the standard 
techniques available in the literature. The ratio of the pt distribution of the electron 
from the decays of heavy flavours produced in heavy ion collisions to the correspond- 
ing (appropriately scaled by the number of collisions) quantities from the pp collisions 
is defined as: R A a(pt) = dN e Au+Au /d 2 p T dy/ [N coU x dN e /d 2 p T dy], which will be 
unity in the absence of re-scattering. The STAR [25J and the PHENIX [26] collab- 
orations have measured the Raa{pt) for non-photonic single electron as a function 
of pt for Au+Au at ^/snn = 200 GeV. The experimental data from both the col- 
laborations show Raa < 1 for pt > 2 GeV indicating substantial interaction of the 
heavy quarks with the plasma particles. The spectra evaluated using the formalism 
described above reproduces the data reasonably well (Fig. [21 right panel). 

Next we discuss the elliptic flow resulting from non-central collisions of nuclei. When 
a heavy quark propagates along the major axis of an ellipsoidal domain of QGP 
(resulting from the non-central collisions) then the number of interactions it en- 
counters or in other words the amount of energy it dissipates or the amount of 
momentum degradation that takes place is different from when it propagates along 
the minor axis. Therefore, the momentum distribution of electrons originating from 
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Fig. 2. Left panel: pt distribution of D mesons. (3 S indicates the value of v r which appears as 
a parameter in the blast wave model. FP stands for the results obtained from the solution 
of FP equation. The experimental data from STAR collaboration |24j is compared with 
the theoretical results. Right panel: variation of Raa with px- The initial temperature and 
thermalization time are taken as 400 MeV and 0.2 fm/c respectively. 
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Fig. 3. Left panel: variation of t>2 with pt for for RHIC energy. Right panel: vt/ 1S 
plotted as a function of pr for LHC energy for 0-10% centrality. The values of Tj and r« 
are taken as 700 MeV and 0.08 fm/c respectively. 

the decays of heavy flavoured hadrons produced from the fragmentation of heavy 
quarks propagating through an anisotropic domain will reflect such anisotropy. The 
degree of momentum anisotropy will depend on both the spatial anisotropy and 
more importantly on the coupling strength of the interactions between the heavy 
quarks and the QGP. The drag and diffusion coefficients depend on the temper- 
ature of the background medium (QGP) which evolves in space and time due to 
expansion. Therefore, the drag and diffusion will also change due to the flow of 
the background. The flow of the background has been treated within the ambit of 
(2+1) dimensional hydrodynamics [27J. The coefficient of elliptic flow, v^ F is defined 
as: V2 F (pr) = (eos(20)) = J dcf)dN/dydp T d(j)\y =0 cos(2(f))/ [J d<pdN/dydpTd(f)\y =0 }. We 
evaluate v^ F in the current formalism [28] and compare the results with experimental 
data [29] (Fig. [3j left panel). The prediction for the elliptic flow of the heavy quarks 
to be measured at Large Hadron Collider (LHC) energy through the semi-leptonic 
decays is depicted in Fig. [3] (right panel). The value of v^ F at LHC is similar to 
that of at RHIC. The prediction for the Raaa at LHC has been displayed in Fig. @] 
separately for D and B mesons. The sensitivity of the results on the equation of 
state (EoS), i.e. on the velocity of sound is also considered. Lowering of c s gives 
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Fig. 4. Variation of Raa with px for LHC for 0-10% centrality. Left panel for charm and 
right panel bottom quarks. The values of Tj and t% are taken as 700 MeV and 0.1 fm/c 
respectively. 

more suppressions as observed in Fig. HI Lower value of velocity of sound, c s makes 
the expansion of the plasma slower enabling the propagating heavy quarks to spend 
more time to interact in the medium and hence lose more energy before exiting from 
the plasma resulting in more suppression. 



4 Summary 



We have shown that the variation of various HBT radii with invariant mass extracted 
from the correlation functions of two lepton pairs can be used to understand the 
evolution of collectivity in ultra-relativistic heavy ion collision. The evolution of the 
radial flow in the system produced in URHIC has been discussed and demonstrated 
that the non-monotonic variation of v r with M can be used as a signal for parton 
to hadron transition. The elliptic flow of the matter probed by the single electrons 
originating from the heavy flavour decays has been studied. The elliptic flow and the 
nuclear suppression factor measured at RHIC have been reproduced and predictions 
for LHC have been given including both the radiative and the collisional processes of 
energy loss in evaluating the effective drag and diffusion coefficients. The sensitivity 
of the results on the EoS has also been studied. 
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